<1177  286  FREOUENCV-DEPENDENT  V-REPRESENTABILITV  IN  DENSITV  '  1/1 

FUNCTIONAL  THEORV(U)  CALIFORNIA  UNIV  SANTA  BARBARA  DEPT 
OF  PHVSICS  D  HEARNS  ET  AL  FEB  87  N88814-84-K-8548 

F/G  12/1 


UNCLASSIFIED 


NL 


AD- A 177  206 


« 


r 


OFFICE  OF  NAVAL  RESEARCH 


Contract  N00014-84-K-0548 
Task  No.  NR372-160 

TECHNICAL  REPORT  NO.  11 


Frequency-Dependent  v-Representability 
in  Density  Functional  Theory 


by 

Daniel  Mearns  &  W.  Kohn 


Prepared  for  publication 
in 

Physical  Review  A  (1987) 


Department  of  Physics 
University  of  California,  Santa  Barbara 
Santa  Barbara,  CA  93106 


DTIC 


Approved  for  Public  Release. 

Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of 
the  United  States  Government. 

Distribution  of  this  Document  is  Unlimited 

February  1987 


87  2  26  048 


OFFICE  OF  NAVAL  RESEARCH 

Contract  N00014-84-K-0548 
Task  No.  NR372-160 

TECHNICAL  REPORT  NO.  11 

Frequency-Dependent  v-Representability 
in  Density  Functional  Theory 

by 

Daniel  Mearns  &  Walter  Kohn 


Department  of  Physics 
University  of  California,  Santa  Barbara 
Santa  Barbara,  CA  93106 


Frequency  dependent  perturbations,  u^(r,u>)e  lti)t,  of  a  system 

in  a  non-degenerate  ground  state  are  studied.  It  is  shown 

that  (l)for  w<oj  .  ,  the  lowest  excitation  frequency  everv 
min 

u,  (r,ai)  causes  a  density  change;  and  (2)  for  to>u)  .  ,  there 
i  min 

are  isolated  frequencies  where  this  is  not  the  case. 


February  1987 


1.1 

V,' 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  rNTiwi  Data  Cnttrtd) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


4.  TITLE  (and  Submit) 


FREQUENCY-DEPENDENT  v-REPRES ENT ABILITY  IN 
DENSITY  FUNCTIONAL  THEORY 


N00014-01 


S.  TYPE  OF  REPORT  A  PERIOD  COVERED 

TECHNICAL  REPORT 
6/21/86  -  6/21/87 


•  .  PERFORMING  ORO.  REPORT  NUMBER 


»■  AU  THORf  t) 


ONTRACT  OR  GRANT  NUMBERS 


DANIEL  MEARNS  &  WALTER  KOHN 


N00014-84-K-0548 


PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

UNIVERSITY  OF  CALIFORNIA 

PHYSICS  DEPARTMENT,  SANTA  BARBARA,  CA  93106 
SEE: CONTRACTS  &  GRANTS -ROOM  3227  CHEADLE  HALL 


II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

OFFICE  OF  NAVAL  RESEARCH 

ELECTRONICS  &  SOLID  STATE  PHYSICS  PROGRAM 

800  N.  QUINCY,  ARLINGTON,  VA  22217 


14.  MONITORING  AGENCY  NAME  A  AOORESSf/t  dllltrtnl  trom  CoruntliHt  Olll et) 

OFFICE  OF  NAVAL  RESEARCH  DETACHMENT 
1030  EAST  GREEN  STREET 
PASADENA,  CA  91106 


program  element,  project,  task 

AREA  S  WORK  UNIT  NUMBERS 


TASK  NO.  NR372-160 


IZ.  REPORT  DATE 


2/6/87 


ts.  NUMBER  OF  PAGES 

3 


IS.  SECURITY  CLASS,  ft/  lb  It  report; 


UNCLASSIFIED 


ISa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (ol  tbit  Rap art) 


"Approved  for  Public  Release:  Distribution  Unlimited" 


17.  DISTRIBUTION  STATEMENT  (•/  the  obotrocl  In  Block  70*  II  different  from  Roper!) 


Reports  Distribution  List  for  0NR  Physics  Division  Office 
Unclassified  Contracts 


IB.  SUPPLEMENTARY  NOTES 


Submitted  to  Physical  Review  A 


If.  ,  KEY  WOSOI:fC«iUnut  on  rower  eo  •  Ido  II  nocoooory  end  Identity  by  block  ntmber) 


-Density  Functional  Theory;  v-Represent ability;  Frequency-Dependent 
Linear  Response; 


20.  ABSTRACT  (Continue  on  rewerme  o Ido  It  nocoooory  end  Identify  by  block  number) 


Frequency  dependent  perturbations,  u^(r,to)e  U)t,  of  a  system  in 
a  non-degenerate  ground  state  are  studied.  It  is  shown  that 
(l)for  ^<UJmin»  the  lowest  excitation  frequency  every  u^(r,io) 

causes  a  density  change;  and  (2)for  to  .  ,  there  are  isolated 

m  in 

frequencies  where  this  is  not  the  case. 


DO  !  j2"*T3  1473  BOITION  OF  1  MOV  «S  IS  OBSOLETE 

S/N  OIOJ-  LF-014-  4601 


SECURITY  CLASSIFICATION  OF  THIS  PAOE  fWNa«  Data  Unttrtd) 


rr 


j 


Fnqaancy>Dap«ad«nt  o-Repre*  mutability 
In  Density  Functional  Theory 


Daniel  Means  and  Walter  Kohn 


Department  of  Phyaiea 
University  of  Celifomie 
Santa  Barbara  CA  93106 


ABSTRACT:  In  density  functional  theory  (DFT)  of  the  ground  state  a  density 
distribution,  rto  (r),  is  called  ^representable  (VR)  if  it  is  the  ground  state  density 
in  some  external  potential.  (It  is  known  that  not  all  “reasonable"  ito(r)  are  VR.)  In 
DFT  of  time-dependent  linear  response  of  a  non-degenerate  ground  state  a  similar 
question  arises:  Is  a  response  density,  VR;  i.e.,  is  it  the  response 

to  some  perturbing  potential  9i(r,u>)e~*w*?  (E.K.U.  Gross  and  W.  Kohn,  Phys. 
Rev.  Lett.  56,  2850  (1985).)  In  the  present  paper  we  show  that  (l),  if  the 
frequency  u  <  u/m,n  (the  lowest  excitation  frequency),  the  answer  is  affirmative; 
and  (2),  if  w  >  u^n,  the  answer  is  not  necessarily  affirmative,  as  demonstrated 
by  counterexamples.  (We  exhibit  “reasonable"  functions  »i(r,u/)  *  which,  at 

isolated  frequencies,  are  not  VR.)  Implications  for  time-dependent  DFT  of  linear 
response  are  discussed. 
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L  INTRODUCTION 


Tie  question  under  what  conditions  &  static  density  distribution,  no  (r) ,  is  unrepresentable 
(VR)1  has  attracted  interest  in  recent  years.2’3’4  The  issue  of  o-representability  also  arises 
for  time-dependant  densities,  n(r,t).  In  particular,  in  connection  with  time-dependent  lin- 
®  response3  one  encounters  the  following  situation:  Let  be  the  non-degenerate  ground 
state  of  a  many  electron  system  with  density  no(r),  in  a  static  external  potential  o0 (r) .  A 
small  perturbing  potential,  t>x(r,  is  known  to  lead  to  a  unique  first  order  density 

response,  rs^(r,fa/)c  ** where  and  are  related  by  the  response  function  x‘ 

rn(r,w)  *  J  x(r,r';w)ei(r/,w)<fr'.  (X) 

The  converse  question  is,  can  a  given  function,  »i(r,w),  be  generated  by  some  function 
wl(r» w)?  This  is  the  o-representability  problem  of  linear  density  response  theory  addnrrr^d 
in  this  paper.  It  may  be  posed  for  interacting  as  well  as  non-interacting  particles. 


IL  NON-INTERACTING  FERMIONS 


For  non-intaracting  fermions,  first  order  time ■  dependent  perturbation  theory  gives  the 
following  result  for  the  response  function: 


=  £(/,•  -  fj) 

ij 


w  -  («j  ”  «i)  +  iS 


(2) 


where  fa  are  the  single  particle  eigenfunctions  of  the  unperturbed  hamiltonian,  et  the 
eigenvalues,  /,•  the  occupation  numbers  (1  or  0)  for  the  ground  state,  and  6  is  a  positive 
infinitesimal 


Except  at  the  resonances,  w  »  —  e^,  x*  is  hermitian  and  real  and  has  a  complete  set 

of  orthonormal  eigenfunctions,  &(r,u/),  and  real  eigenvalues,  A^(u/): 

J  Xs(r,r';w)ft(r/Ia;)dr/  =  A*(w)&(r,w).  (3) 

S’,  for  some  frequency  Q,  one  of  tho  eigenvalues,  say  Aj,  vanishes  then  the  perturbing 
potential 

0l(r,O)  3^(r,i5)  (4) 

{n  «  1)  has  vanishing  first  order  density  response: 

*l(r,fi>)  =  J  Xj(r,r';<S)01(r',O)dr' 

i 

*  mA *(6>)$(r,u>)  =  0.  (5) 

Also,  clearly,  the  density  change 


r»i(r,w)  =  nq{r,Q)  (6) 

is  not  induced  by  any  linear  combination  of  the  complete  set  &(r ,£),  that  is,  it  is  not  VR 
as  a  linear  response  density. 
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Conversely,  if  ail  Xg(u)  0  any  density  that  can  be  expressed  by  a  series  of  the  complete 

set  of  eigenfunctions, 

*1  *  £  »l,*(w)fc(*’*«).  (7) 

l 

is  VR  by  the  potential 

t>l(r,w)  =  JlA<(a;)nu(w)a(r,a/),  (8) 

l 

provided  that  the  latter  series  converges. 

For  all  frequencies  there  exists  one  vanishing  eigenvalue,  corresponding  to  the  trivial 
perturbation  vi  =  const.  (The  corresponding  «i  *  const,  is  trivially  not  VR.)  We  shall 
show  that  for  any  system  at  frequencies  smaller  than  the  first  resonance  no  non-trivial  van¬ 
ishing  eigenvalues  exist,  and  hence  all  densities  which  have  a  sufficiently  rapidly  convergent 
expansion  in  the  functions  &  are  VR. 

Restricting  attention  to  non-resonant  frequencies  and  choosing  the  eigenfunctions  <£,• 
real,  Eq.  (2)  may  be  written  as 

X.frr'ju.)  -  £  *.(r)*#(r)^(r')*«(r')t  (9) 

aj  W  €0a 

where  the  indices  a  and  0  denote  the  occupied  and  unoccupied  levels,  respectively.  The 


response  is 


where 


»i(*,.w)  =  £  2  t0aj~  vfaM*)MT)' 

a,0  W  *£<* 


V0a  5  J  40(r)vi(r,u)4a(r)dr.  (11) 

The  matrix  elements  Vpa  cannot  all  vanish,  otherwise  the  determinantal  wavefunction 
perturbed  by  would  be  identical  to  the  unperturbed  ground  state.  The  integral 

f  »»i  (r,  w)vi  (r,  u)dr  *  vja  (12) 


KIPS 


is  then  negative-definite  for  frequencies  below  the  first  resonance.  Therefore,  n\  can  not 
vanish  identically,  and  Xs  can  have  no  vanishing  eigenvalue  in  this  range. 

Let  us  note  in  passing  that,  for  the  special  case  of  a  single  particle,  the  density  change 

is 

»l(r,w)  =  ^i(r)  J2  (13) 

Owing  to  the  linear  independence  of  the  eigenfunctions,  can  not  vanish  identically 

at  any  frequency. 

We  shall  now  present  two  examples  of  systems  which,  at  isolated  frequencies  above  the 
first  resonance,  have  non-VR  response  densities. 


A.  One-dimensional  Ring,  00=  0 

For  a  one-dimensional  ring,  0  <  x  <  2r,  the  non-interacting  eigenfunctions  are  plane 
waves: 


*<(*)  =.  (2t)'V2  «“**. 


(14) 


In  the  rrwmnrm  gauge,  Jk*  is  an  integer.  The  ground  state  with  one  particle  is  of  no 
interest  (always  VR).  The  ground  state  with  two  particles  is  two-fold  degenerate  and  hence 
inadmissable  for  our  purposes.  However,  by  choosing  an  appropriate  constant  gauge,  the 
allowed  kt  are  all  shifted  by  \  so  that  kt  =  [21  +  l)/2,  t  *  0,  ±1,  ±2, ....  In  this  gauge,  the 
two-particle  ground  state  is  non-degenerate.  Eq.(2)  takes  the  form 


!  ei(k9-ka)x  e-i(k9-ka)x/  e~i(k9-ka)x  ei(k9-ka)x' 

—  2jt  ^  ^  w  -  60a  U  +  (0a 

<*,P 


(15) 


where 


(16) 


5 


The  eigenfunctions  of  Xs  are  plane  waves,  independent  of  u/,  so  that 


where 


and 


Xs  =  J2 
#0 


aw  =  (2*rl/2  *“*, 


5  w  -  (*1  -  *2)  5  ‘•'+(*?-*2) 


(17) 


(18) 


-  L 


V  w-f^  +  Mk.)  u+(*>-2tt„)' 


(19) 

with  iV  the  number  of  particles.  For  example,  for  JV  *  2,  *  1  and  A:a  =  so  that 

A<<")  “  +  +  ZT-(ei  -if’  1 4 1  >  !•  (“) 

Hence  each  \g,  J£J  >  1,  has  two  poles,  with  a  zero  lying  between  them  at  <3  *  £(£2  -  l)1/2. 


B.  One-dimensional  Box,  v0=  0 

For  a  one-dimensional  box,  0  <  x  <  r,  the  non-interacting  eigenfunctions  are  standing 
waves:  l 

Mx)  =  (2/w)1/2sin£x,  (21) 

where  £  =  1,2, 3,...,  and  so 


X$(x,x';u)  =  —  Y'  -0  6/?a — s-  sin  ax  sin/Jx  sin/3x/  sinax', 


where 


-  a2. 
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(22) 


(23) 


With  the  expansion 


(24) 


oo 

t>i(x,u/)  =  a((u)coelx, 
t=l 

the  matrix  elements,  Eq.  (11),  take  the  form 

y*  =  2  (a^— a  ”  a/?+a)» 


hence 


(°0-<*(w)  -  <*3+a(w))  [co*(^  -  a)*  -  c<*(0  +  *)*]• 


The  eigenvalues  and  eigenfunctions  of  x*  ere  given  by 


i 

t 


(25) 


(26) 


i>i(z,u)  m  A(u)si(z,u). 


(27) 


For  two  particles  this  leads  to  the  set  of  equations: 

Aai  “  ;  <•>  -  *»>• 

a°2  =  \  {^rzgj  (“2  - •*) +  (•»  -  •*)} 

(“5 -<■»)  +  (a3 -“7)} 

Aa4  °  z  (s,2  -  242  («6  -  «6)  +  ^  ”3^2  (<M  -  ««)  («2  -  <M) } 

1  r  (*  +  i)2-i  ,  x  .  (£  +  2)2  —  4 

Aa* "  ;  u-((/+i)»-i|2 {at " aM) + ^-[(f+V-ii2  ( ' ' at+4) 

i  [(1  l)2t)7  lip  (^-2  -  0  ~U2  i  (d  tf-  4|2  -a;>};  '-5-  (28) 


Since  the  even  and  odd  Fourier  components  are  not  coupled  the  eigenfunctions  have 
definite  parity.  Numerical  solutions  of  the  equations  corresponding  to  even  eigenfunctions 
have  been  carried  out  for  frequencies  0  <  u>  <  50.  The  method  consists  of  taking  a  finite 
series  for  Eq.  (24),  so  that  the  solution  of  Eqs.  (28)  is  reduced  to  diagonaiization  of  a 
finite  matrix.  Since  for  any  fixed  u /  the  series  of  Eq.  (22)  is  uniformly  convergent  in  the 
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variables  x  and  x/,  Xs  can  be  approximated  with  arbitrary  accuracy  by  a  sufficiently  large 
matrix,  for  any  fixed  frequency  range. 

A  plot  of  the  eigenvalue*  for  a  30-dimensional  matrix.  Fig.  1,  exhibit*  two  eigenvalue* 
passing  through  zero,  above  the  resonances  at  12  and  32.  Fig.  2  is  an  expanded  view  of  the 
first  zero  in  which  the  repulsion  between  the  eigenvalues,  indicating  a  mixing  between  the 
eigenfunctions  in  regions  of  near-degeneracy,  is  more  pronounced.  Despite  this  mixing  the 
eigenfunction  corresponding  to  the  vanishing  eigenvalue  tends  to  a  limit  as  the  eigenvalue 
approaches  zero  from  above,  as  shown  in  Fig.  3  for  the  first  zero.  Fig.  4  is  a  plot  of  the 
eigenfunction  at  a  frequency  for  which  the  eigenvalue  is  very  small  and  positive.  It  is  given 
by  a  rapidly  converging  Fourier  series  and,  accordingly,  is  smooth  in  appearance. 


EEL  INTERACTING  FERMIONS 


The  response  function  for  N  interacting  fermions  is 

2£  r 

x(r,p';w)  =  ^  2  _kl2  /  ^<’(r»r2i‘">r.v)^t(r,r2,---,riv)<fr2,-'<il‘v 

4=1  w  ^ko  J 

x  f  ♦k(r/.r/2»--*.r,^),iro(r/)r/2,---.r/iv)^,/2--T/Jv,  (29) 

where  2?t  and  'F,-  (0  <  i  <  oo)  are  the  eigenvalues  and  normalized  eigenfunctions  of  the 

JV-particle  hamiltonian,  and  =  Eh  —  Eq.  The  relations  analogous  to  Eqs.  (10)  and 
(12)  for  the  non-interacting  case  are 

r»i(r,u;)  -  vko  /''iro(r,r2,---,r^)^rfc(r,r2,---,r^v)dr2---driv  (30) 

4=1  *ko  J 

and 

/  »i(».«)n  (».«)*  =  £  Vfc2°’  (31) 

4=1  "c,4o 

where 

Vko  =  M  J  ^4(rl>r2»*  *  *rAr)t,l(ri)^o(rl»r2> ' '  •  >r^)dridT2  •  •  •  dr^y.  (32) 

Eq.  (31),  like  Eq.  (12),  is  negative-definite  for  frequencies  below  the  first  resonance. 
Therefore,  exactly  as  shown  in  Sec.  I  for  Xai  it  follows  that  x  can  have  no  vanishing 
eigenvalues  in  this  range. 


IV.  CONCLUDING  REMARKS 


In  their  paper  on  density  functional  theory  of  linear  response.  Gross  and  Kohn5  pre¬ 
supposed  that  the  physical  density  rio(r)  +  (r,  t)  was  “non-interacting  VR"  (VR-N), 

that  is,  can  be  reproduced  by  a  system  of  non-interacting  particles  in  an  external  potential 
v0(r)  +  vi  (r,  t).  We  have  shown  in  this  paper  that  this  will  be  the  case  if  by  itself  is  VR- 
N  and  the  frequency  of  vj  is  less  than  the  smallest  resonance.  However,  if  the  frequency 
is  higher,  our  examples  show  that  caution  is  in  order.  If  our  examples  are  representative, 
then  in  general  we  expect  that  there  will  be  isolated  frequencies,  <2,  at  which  most  density 
changes  are  not  VR-N.  The  exceptions  are  those  special  functions  which  are  orthogonal  to 
the  functions  fjr(r,u))t  corresponding  to  vanishing  eigenvalues  of  x,. 
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FIGURE  CAPTIONS 


Fig.  1.  Eigenvalue*  of  x<  for  *  30-dimensional  matrix.  The  dashed  lines  mark  the 
locations  of  the  resonances  at  w  =  8,  12,  24,  32,  and  43. 

Fig.  2.  Eigenvalues  of  x*  in  a  frequency  range  containing  a  aero. 

Fig.  3.  Magnitudes  of  the  first  10  Fourier  amplitudes  of  an  eigenfunction  with  eigen¬ 
value  t^rwiitig  to  aero,  weighted  so  that  (r/2)1//2^|a^|2  =  1.  The  dashed  line  is  the 
eigenvalue  curve  in  the  same  frequency  range  (vertical  scale  not  shown). 

Fig.  4.  Normalised  eigenfunction  for  X(u)  *  1.2386xl0~5,  w  =*  12.658. 


